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Abstract. In non-relativistic quantum mechanics, singular potentials in problems with 
spherical symmetry lead to a Schrodinger equation for stationary states with non-Fuchsian 
singularities both as r — > and as r — > oo. In the sixties, an analytic approach was 
developed for the investigation of scattering from such potentials, with emphasis on the 
polydromy of the wave function in the r-variable. The present paper extends those early 
results to an arbitrary number of spatial dimensions. The Hill-type equation which leads, 
in principle, to the evaluation of the polydromy parameter, is obtained from the Hill 
equation for a two-dimensional problem by means of a simple change of variables. The 
asymptotic forms of the wave function as r — > and as r — > oo are also derived. The 
Darboux technique of intertwining operators is then applied to obtain an algorithm that 
makes it possible to solve the Schrodinger equation with a singular potential containing 
many negative powers of r, if the exact solution with even just one term is already known. 
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1. Introduction 

One of the longstanding problems of non-relativistic quantum mechanics is the investi- 
gation of scattering from singular potentials, with efforts by many authors over several 
decades (see [1-15] and references therein). The main motivations can be described as 
follows [13]. 

(i) Repulsive singular potentials make it possible to obtain a fairly accurate description of 
the short-range part of the nucleon-nucleon interaction. 

(ii) The (p,p) and (p,n) processes can be interpreted in terms of complex potentials r~ n , 
with n > 2. 

(hi) Repulsive singular potentials reproduce also the interactions of nucleons with K- 
mesons, and a — a scattering processes. 

(iv) The Lennard-Jones potential, proportional to r~ 12 , can be used to study interactions 
among the overlapping electron clouds of non-polar molecules. 

(v) At a field-theoretical level, it appears quite remarkable that non-renormalizable field 
theories give rise to effective potentials in the Bethe-Salpeter equation which are singular 
[1, 2], whereas super-renormalizable and renormalizable field theories give rise to regular 
or transition-type effective potentials, respectively. There was therefore the hope that any 
new insight gained into the analysis of non-relativistic potential scattering in the singular 
case, could be eventually used to obtain a better understanding of quantum field theories 
for which perturbative renormalization fails (cf section 6). 

(vi) In particular, one might then hope to be able to "map" the analysis of quantum gravity 
based on the Einstein-Hilbert action (plus boundary terms), which is well known to be 
incompatible with the requirement of perturbative renormalizability [16], into a scattering 
problem in the singular case, for which the Schrodinger equation for stationary states: 




(1.1) 
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has non-Fuchsian singularities (see appendix) both as r — > and as r — > oo. With our 
notation, g is the number of spatial dimensions, and 1(1 + q — 2) is obtained by studying 
the action of the Laplace-Beltrami operator on wave functions belonging to the tensor 
product [17] 

L 2 (K + ,r q - 1 dr)®L 2 (S q - 1 ,dn). 

Moreover, with a standard notation, one has (of course, the energy E is positive in a 
scattering problem) 

k^ 2 -f (1.2) 

V(r) ee ^U{r) (1.3) 

with £7 (r) the potential term in the original form of the Schrodinger equation. From now 
on, it is V(r) which will be referred to as the potential, following the convention in the 
literature. 

Among the analytic results obtained so far in the investigation of potential scattering 
in the singular case, we find it appropriate to mention what follows. 

(i) A constructive determination of the S'-matrix, based on the polydromy properties of 
the wave function (see appendix) and on the Hill equation for the polydromy parameter 
[6, 8]. 

(ii) Perturbative technique for the potential V(r) = g 2 r~ 4 in three dimensions, by re- 
expressing the radial Schrodinger equation as a modified Mathieu equation [10], with eval- 
uation of S'-matrix and Regge poles. 

(iii) Generalized variable-phase approach, leading to a JWKB phase-shift formula [13]. 

(iv) Generalization of the JWKB method to arbitrary order, with rigorous error bounds 
[14]. 

In our paper, sections 2 and 3 apply the method of [6, 8] to the Schrodinger equation 
for stationary states in three or more spatial dimensions, proving that a simple but deep 
relation exists between the corresponding Hill equations in two and three or more spatial 



dimensions. Section 4 presents, for completeness, the JWKB analysis of the wave function, 
jointly with its limiting behaviour as r — > and as r — > oo. Section 5 studies the application 
of the intertwining operator technique to singular potential scattering. Results and open 
problems are described in section 6. 

2. Schrodinger equation for stationary states 

Following the remarks in the introduction, we first study the Schrodinger equation for 
stationary states in three spatial dimensions in a central potential: 



d? | 2 d 1(1 + 1) | fc2 
dr 2 r dr r 2 



tP(r) = V(r)ifj(r). (2.1) 



What is crucial is the polydromy of the wave function in the r variable. Indeed, if the 
potential V(r) is a single-valued function of r, one can find two independent solutions 

Mr) = r^xi(r) (2.2) 

Mr) = r" 7 X 2(r) (2.3) 

where xi an d Xi are single-valued functions of r, and 7 is a parameter which can be 
determined from a transcendental equation (see below). The general solution of Eq. (2.1) 
is therefore of the form 

<0(r) = oil ipi(r) + a 2 fair)- (2.4) 

Remarkably, one can compute directly Xi( r ) an d an d study their behaviour as r — > 

and as r — > 00 [6, 8]. For this purpose, the following Laurent expansions are used (the 
subscript for x is omitted for simplicity): 

00 

W(r) = r 2 [V(r) - k 2 ] = w n r n r e]0,oo[ (2.5) 

n= — 00 

00 

X(r)= ^ re ]°'°°[- ( 2 - 6 ) 
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These expansions hold because V(r) is assumed to be an analytic function in the complex-r 
plane, with singularities only at infinity and at the origin [6, 8]. The Laurent series (2.5) 
and (2.6) are now inserted into Eq. (2.1), which is equivalent to the differential equation 

(cf [6, 8]) 

r 2 |§ + 2( 7 + l)r^f+( 7 (7 + 1) - *(* + l))x = r 2 [v(r) - k 2 ] X . (2.7) 
One thus finds the following infinite system of equations for the coefficients (cf [8]): 

oo 

(n + 7)(ra + 7 + 1) - A c n = ^ 



m= — oo 



where 

X 2 = l(l + l) + w (2.9) 

u n = w n - w S nfi - (2-10) 

To solve the system (2.8) one first writes an equivalent system for which the determinant 
of the matrix of coefficients is well defined. Such a new system is obtained from (2.8) by 

2 

dividing the n-th equation by (n + 7 )(ra + 7 + 1) — A . The resulting matrix of coefficients 
has elements 

tt _S Un-m (0 M\ 

^n^m — u n,m 2 

[(ra + 7)(n + 7+ 1) - A ] 
where det(H) exists since the double series 



E 



[(n + 7 )(n + 7 + l) - A 2 ] 



converges for all values of 7 which do not correspond to zeros of the denominator. At this 
stage one can appreciate the substantial difference between regular and singular potentials. 
In the former case, u n is non-vanishing only for positive n. In the singular case, however, 
the presence of negative powers in the Laurent series (2.5) gives 7 as the solution of a 



5 



transcendental equation, i.e. (the vanishing of det(H) being necessary and sufficient to 
find non-trivial solutions of the system (2.8)) 



F( 7 ) = det(H) = 0. 



(2.12) 



3. Equation for the 7 parameter 



To evaluate ^(7), we point out that, on defining 



7 = 7 + 



(3.1) 



A 2 = A" + 7 =(/ + -!- ) +w 



-2 1 



one finds 



—2 



(n + 7)(n + 7 + 1) - A~ = (n + 7 ) 2 - A 2 . 



(3.2) 



(3.3) 



This simple but fundamental property makes it possible to perform the three-dimensional 
analysis by relying entirely on the investigation in two spatial dimensions, because H n ^ m 
now reads 

Un—m 



^ n ' m " 5n ' m j(n + 7) 2 -A2] 
and hence, from the work in [6, 8], one knows that 



(3.4) 



F(t) = 1+ F(0) - 1 



cot 7r(7 + A) — cot 7r(7 — A) 



2 cot TV A 



(3.5) 



where F is an even periodic function of 7, with unit period [8]. The equation 



F(7) = 



(3.6) 



is, as we said in section 2, a transcendental equation. If a root, say x, is known, at least 
approximately, one can then evaluate the desired 7 parameter from the definition (3.1) as 



7 = x — 



1 



(3.7) 



The ground is now ready for understanding the key features of singular potential 
scattering in an arbitrary number of spatial dimensions. For this purpose, we remark that, 
upon setting ip(r) = r 7 x(r), Eq. (1.1) leads to the following second-order equation for x 
(cf(2.7)): 

o d 2 _ N d 



dr 2 



+ (2j + q-l)r 



dr 



+ ( 7 2 + (q- 2) 7 -1(1 + q- 2))] X (r) = W(r) X (r). 



(3.8) 



Thus, on defining (cf (3.1)) 



7 = 7+2(9- 2 ) 



(3.9) 



one can re-express Eq. (3.8) in the form 

o d 2 \ d 



dr 2 



(27+1)' 



+ (f- ^-2) 2 -/(/ + g-2)) 



x (r) = W(r)x(r). 



(3.10) 



At this stage, the Laurent expansions (2.5) and (2.6) lead to an infinite system of equations 
for the coefficients c n in the form (cf (2.8)) 



(n + 7) 2 -A s 



(3.11) 



m= — 00 



where we have defined (cf (3.2)) 



A 2 /(/ • q 2) • X -\n 2) 2 +w =(l+ ^{q - 2) ) + ,r {] 



(3.12) 



whereas the notation (2.10) remains unchanged. Thus, one can always perform the analysis 
in terms of the infinite matrix (3.4), provided that one defines 7 and A 2 as in (3.9) and 
(3.12), respectively. The resulting Hill-type equation which leads, in principle, to the 
evaluation of the fractional part of the polydromy parameter 7, involves an even periodic 
function of 7 + \(q — 2). 
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4. Asymptotic form of the solutions 

Since one might be eventually interested in the S'-matrix, it is quite important to study 
the limiting behaviour of stationary states as r — > and as r — > oo. In the former case, 
one can perform a JWKB analysis of Eq. (1.1), setting therein 

if>(r) = A(r)e^ 5(r) . (4.1) 
This leads to the equation (the prime denoting differentiation with respect to r) 
(2m(E - U(r)) - S' 2 ^A + ih (2A S' + AS" + ^-y^AS'^J 



+ v[A" + ^A>- l « + y 2 h 



= 0. (4.2) 



If, in a first approximation, the term on the second line of Eq. (4.2) is neglected, one finds 
the equations 

S' 2 = 2m(E-U(r)) (4.3) 



±L(A 2 S') + ^—^-A 2 S' = (4.4) 
dr r 



which imply 



S' = ±^2m(E-U(r)) (4.5) 

A 2 S' = constant x r"^" 1 ) (4.6) 

and hence, for some constant 

_ i 

A(r) = Pr- i3 T 1 (2m(E-U(r)j) ~\ (4.7) 

To second order in h, one has to consider the second line of Eq. (4.2). On taking the 
prefactor A(r) in the form (4.7), one has then to evaluate the phase S(r) from the equation 

S'(r) = ±^2m(E-U(r)) + h 2 f lq (A(r)) (4.8) 
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where 



fl,(A(r)) = ^1 + (l ~ ^ A ' + 



1 



A 



(q 2 -4q + 3) + l(l + q-2) 



+ Vlu" (2m(E - U(r)j) * + ^m 2 U' 2 (2m(E - U(r)j) ' 



Only an approximate calculation of the square root in Eq. (4.8) is possible, if 



P 2m(E-U(r)) 



(4.9) 



by expanding yjl + p in powers of p, but this does not improve substantially the under- 
standing of the behaviour of the wave function as r — > for a fixed value of k (see below), 
and hence we do not present further calculations along these lines. One should bear in 
mind, however, that the JWKB expansion has an asymptotic nature, and rigorous error 
bounds can be obtained [14] . 

In particular, in the physically more relevant case of three spatial dimensions, Eq. 
(4.7) leads to (see (1.2) and (1.3)) 



_ i 

A(r)=pr- 1 (k 2 -V(r)^j \ 



(4.10) 



When r — > 0, V(r) is much larger than k 2 for a fixed value of k, and hence the JWKB 
ansatz (4.1) leads to (hereafter U = 1) 



~ Bjjj r _1 (y(r))"4 exp / \JV{y) dy 

J r 



(4.11) 



for some parameters Bjjj depending on 7 and / (cf [6, 8, 18]). Of course, the JWKB 
solution for all values of k which results from (4.5) and (4.10) is, instead, 



ipi,n ~ Pi,nr (k - V(r)) 4 exp 
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V^ 2 - V(y) dy 



ro 



(4.12) 



Moreover, as r — > oo, one has the familiar asymptotic behaviour 

i>i,n ~ A+ij r" 1 exp jz J/cr - Z^] j 

+ AJjj r _1 exp{-i [fcr-Z^]} . (4.13) 

The S'-matrix is given by the formula [6, 8] 

(a+B h - A+BA 

S = j ( (4.14) 

[AjBn - A-Bjj 

where the A and B parameters are the ones occurring in the asymptotic expansions (4.11) 
and (4.13), and can be obtained by means of the saddle-point method [6, 8]. 

5. Intertwining operators for singular potentials 

Since exact solutions of singular scattering problems in terms of special functions are known 
in a few cases only, it appears quite important to look for a technique that makes it possible 
to generate solutions for complicated problems, relying on what is known in simpler cases. 
For this purpose, we here consider the Darboux method of intertwining operators [19-22]. 

The aim of the Darboux method is to generate families of isospectral Hamiltonians. 
It relies on a theorem which, in modern language, can be stated as follows [23]. Let ip be 
the general solution of the Schrodinger equation 



Htj}{x) = 



d 



2 



+ V(x) 



4>(x) = Eip(x). (5.1) 



dx 2 

If (p is a particular solution of (5.1) corresponding to an energy eigenvalue e ^ E, then 

j = l (5 .2) 



<p \ dx dx 
is the general solution of the Schrodinger equation 



H^{x) = E^{x) (5.3) 
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where 



H = 



dx 2 



+ V(x) 



(5.4) 



V(x) = V(x) -2— log ip(x). 



(5.5) 



In other words, if two Hamiltonian operators, say Ha and Hb, are given, one looks 
for a differential operator, say D, such that [24] 



H B D = D H A . 



(5.6) 



It is then possible to relate the eigenfunctions of Ha and Hb by using the action of D (see 
below). Here, we focus on one-dimensional problems, with 



Ha = Hn 



Hb = Hi = 



d 



d 2 



dx 2 



dx 2 



+ V (x) 



+ V 1 (x) 



D = Tx +G ^ 



(5.7) 



(5.8) 



(5.9) 



where Vo and V\ are the "potential" functions, and G is another function, whose form is 
determined by imposing the condition (5.6). This reads, explicitly, 



d 2 



dx 2 



+ V 1 



dx ^ / \ dx ^ 



d 2 



dx 2 



+ Vo)f 



(5.10) 



for all functions / which are at least of class C 3 . On imposing Eq. (5.10), one finds exact 

cancellation of the terms — ^! and — since they occur on both sides with the same 

sign. Hence one deals with the equation 



( - 2G' + V 1 - Vo) ^+( - G" - V ' + (V! - V )g) 



f = 



which implies 



2G' = V 1 -V Q 



(5.11) 



(5.12) 
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—G" - Vq + (Vi - Vq)G = (5.13a) 

by virtue of the arbitrariness of /. It is now possible to use Eq. (5.12) to express Eq. 
(5.13a) in the form 

+ (5,3.) 

which is solved by 

G 2 - G' = Vq + C (5.14) 

for some constant C. Equation (5.14) is known as the Riccati equation. Its non-linear 
nature makes it desirable to develop an algorithm to relate it, instead, to the solution of a 
linear problem. This is indeed achieved by considering the function ip such that 

G=~logip. (5.15) 

The Eqs. (5.12) and (5.15) are, of course, in complete agreement with the result (5.5), 
with V replaced by Vi, and V replaced by Vq. One then finds, by virtue of (5.14) and 
(5.15), that cp obeys the linear second-order equation 

H ip = -Ccp. (5.16) 

This is a simple but deep result: one first has to find the eigenfunctions of Hq, say ip, 
belonging to the eigenvalue — C. Once this is achieved, the desired function G is obtained 
from (5.15), and hence the intertwining operator is 

D = ±-*. (5.17) 

ax <p 

In the applications, it is also convenient to use Eq. (5.12) to express Eq. (5.14) in the form 

G 2 + G' = V 1 + C. (5.18) 

If one studies Eq. (1.1) or, in particular, Eq. (2.1), the standard definition in three 
dimensions 

M r ) = (5.19) 
r 
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leads to a second-order differential operator acting on y which is of the form (5.7) or (5.8). 
However, the choice of a suitable intertwining operator, aimed at relating operators Ha 
and Hb whose potential terms differ in a somehow substantial way, is a non-trivial task. 
For example, if one considers 

+ f (5-20) 

Eq. (5.18) may be then satisfied by (cf (5.9)) 

G(r) = A (5.21) 

provided that A = k 2 , B = —2k and C = 0. However, the resulting potential Vb(r) is 
found to be, from Eq. (5.12), 

k 2 9 k 

Mr) = =r + f (5.22) 

so that the intertwining operator ends up by relating operators Ha and Hb whose potential 
terms have precisely the same functional form. A scheme of broader validity, however, is 
obtained by looking for V\{r) and G{r) in the form of (Laurent) series, i.e. 

V l (r)= a nr n (5-23) 

n= — oo 
oo 

G(r) = b P rP - ( 5 - 24 ) 

p= — oo 

The insertion of (5.23) and (5.24) into Eq. (5.18) leads to the infinite system 

oo 

(n + l)b n+ i + ^ b p b n - p = a n + C5 nj0 (5.25) 



p= — oo 



that should be solved, in principle, for b n , for all n. One then finds, from Eq. (5.12), a 
Laurent series for V as well, i.e. 



oo 



Vo(r)= ^ ( 5 - 26 ) 



n= — oo 
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where 

oo 

fn = -{n + l)b n+1 + ^ b p b n- P ~ C5 nfi - (5.27) 



p= — oo 

,2^.-4 



For example, if one takes V\(r) = g r , one has 

a n = g 2 S n -4 (5.28) 
and hence one deals with the infinite system 

oo 

(n + l)b n+1 + b P b 

n—p g 3n,— 4 + C5 ni0 . (5.29) 

p= — oo 

This is a non-linear algebraic system for which it does not seem possible to obtain a 

solution such that only a few b p coefficients are non- vanishing. For example, if one tries to 
get bp = unless p = —3, —2, —1, one finds, on setting n = —6, —5, —4, —3, —2, in (5.29) 
the six equations 

(b- 3 ) 2 = (5.30) 

26_ 3 6_2 = (5.31) 

-36_ 3 + 26_ 3 6_i + (6_ 2 ) 2 = g 2 (5.32) 

2(6_i - 1)6-2 = (5.33) 

(6_i - l)6_i = (5.34) 

= C (5.35) 

whereas n = — 1 leads to a trivial identity. Now equations (5.30) and (5.31) imply that 
6_3 = 6_2 = 0, and hence g 2 = from (5.32), which is incompatible with our assumptions. 
The remaining equations (5.33)-(5.35) allow for 6_i = 1, further to 6_i = 0, but with 
C = 0. 

However, the implications remain of high interest: to find non-trivial solutions with 
g 2 7^ and C ^ one needs a large number of b p coefficients (maybe infinitely many), 
including those with p > 0. This still means that one has the opportunity to solve the 
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Schrodinger equation with a complicated singular potential, starting from what one knows 
when the potential equals g 2 r~ 4 [8]. For this purpose, on denoting again by 9? the eigen- 
function of Hq belonging to the eigenvalue — C, and by x = Dtp the eigenfunction of Hi 
belonging to the same eigenvalue, we notice that the desired 9? can be written in the form 



where K(r,r') denotes the Green kernel of the intertwining operator D = -f - + G(r). We 
need such an integral formula because we have chosen, in our particular example, the 
form of the potential term V\(r) in the Hamiltonian H±, for which the scattering states 
are already known in the literature [8]. The unknown are instead the scattering states 
resulting from the Hamiltonian operator with potential term equal to Vb(r) (see (5.26)). 

6. Concluding remarks 

Our paper has studied some aspects of scattering from singular potentials in quantum 
mechanics. Its contributions are as follows. 

(i) The technique of Fubini and Stroffolini [6, 8] , with emphasis on the polydromy properties 
of the wave function, has been applied to an arbitrary number of spatial dimensions, 
say q, when the potential admits a Laurent series expansion. The equation obeyed by 
the polydromy parameter, 7, involves a function which is an even periodic function of 
7+^((/ — 2). Interestingly, one can rely entirely on the analysis performed in [6, 8], provided 
that one considers the parameters defined in (3.9) and (3.12) (strictly, the authors of [6, 8] 
start from three dimensions, but use a transformation [9] leading to an equation formally 
analogous to the radial part of the stationary Schrodinger equation in two dimensions). 

(ii) The Darboux technique of intertwining operators has been applied to relate the singular 
potential terms in the Schrodinger equation for stationary states. The algorithm resulting 
from Eqs. (5.23)-(5.27) leads, in particular, to the non-linear algebraic system (5.29) if 
the potential V\ is taken to be g 2 r~ 4 . 




(5.36) 
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Ultimately, one might want to use these properties to study quantum field theories 
which are not perturbatively renormalizable, according to the original motivations for this 
research field [1, 2, 9]. For this purpose, it seems crucial, to us, to consider the quantum 
gravity problem, focusing (at least) on the following questions. 

(i) What is the counterpart, in quantum gravity, of the Bethe-Salpeter equation containing 
effective potentials of the singular type? As is well known, this equation arises in the 
course of studying the quantum theory of relativistic bound states, and unfortunately a 
simple extension of the Schrodinger equation is not available [25]. Even on neglecting 
curvature effects due to gravitational fields, one then faces retardation effects which lead 
to an extra relative time variable in the problem [25]. An alternative description uses 
a mediating field, whose quantum properties, however, cannot be ignored [25]. When a 
quantum theory of gravity is considered in a space-time approach [26], one may expect 
of using the (formal) theory of the effective action, with a corresponding set of integro- 
differential equations. These should be solved, in principle, by using the functional calculus. 
But even if one were able to achieve so much, and hence derive an effective potential which is 
a gravitational counterpart of the potential normally used to reduce the number of degrees 
of freedom of relativistic bound-state problems, the problem would remain of giving a 
proper interpretation of such potentials, since they seriously affect the exact theory and 
may introduce fictitious singularities (see page 493 of [25]). 

(ii) What kind of results can be then "imported" on mapping quantum gravity into a 
scattering problem from singular potentials (e.g. the asymptotic behaviour of the phase 
shift [9, 13], the exact or approximate solutions derived with some particular choices of 
singular potentials [1-8, 10, 14], the existence theorem for the wave operators [12])? 

(iii) How fundamental is the Darboux method of intertwining operators [19-22] proposed 
in our section 5? Since also the variable phase approach to potential scattering relies on 
an equation of Riccati type [7] , such a question appears to be non-trivial. 

The above issues seem to point out that new perspectives are in sight in the analysis 
of potential scattering for a wide class of singular potentials, with possible implications for 
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a longstanding problem, i.e. the key features of a quantum theory of the gravitational field 
(see [16] and references therein). Hence we hope that the present work, although devoted 
to some technical issues, may lead to a thorough investigation of quantum gravity from a 
point of view well grounded in the general framework of modern high energy physics (cf 
[27]). 

Acknowledgments 

The author is much indebted to Giuseppe Marmo for teaching him all he knows about the 
Darboux method. The present paper is dedicated to the memory of Roberto Stroffolini. 

Appendix 

To be self-contained, let us describe what is meant by Fuchsian singularities of second-order 
differential equations. A theorem due to the German mathematician Immanuel Lazarus 
Fuchs states that a necessary and sufficient condition for the linear equation 



d 2 , N d 



+ Pi(x)— +p 2 (x) 



dx 2 dx 



y(x) = (Al) 



to admit a fundamental system of integrals, say yi(x) and y2(x), which, in the neighbour- 
hood of the singular point xq, can be expressed as [with ^2,^2,^ analytic functions in the 
neighbourhood of xo, for some constants r±, r 2 , fi, A : {r\ — r 2 ) ^ Z] 

yi(x) = (x - x ) Tl (pi(x) (A2a) 

y 2 (x) = (x-x ) r2 (p2(x) (A2b) 

or 

yt(x) = (x - x ) ri (pi(x) (A3a) 

V2{x) = yi(x) A\og(x - x ) + ip(x) (A3b) 

is that pi and p 2 should have poles of order < 1 and < 2, respectively, at the singular 
point xq. One then says that xq is a Fuchsian singularity for equation (Al). 
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To study the point at infinity, one defines 



X 



(A4) 



which leads to the equation (cf (Al)) 



+ 



d 1 fl 



y(0 = o. 



(A5) 



In the analysis of equation (A5) as £ — > 0, which corresponds to the point at infinity of 

(Al), one can thus use again the Fuchs theorem, which implies that p\ (^j and pi (^jj 

should have zeros of degree > 1 and > 2, respectively, at £ = 0. If this condition is fulfilled, 
equation (Al) is said to be Fuchsian at infinity. 

When all singular points are Fuchsian, the corresponding differential equation is said 
to be totally Fuchsian. The non-Fuchsian singularities are, by contrast, singular points of 
(Al) for which the above conditions on poles and zeros of the functions p\ and P2 are not 
fulfilled. 

In our paper, the word polydromy refers to the well known property of some functions of 
complex variable of being multi- valued functions of the independent variable. For example, 
if z is complex, its logarithm is given by the formula 



log(z) = log \ z\ + % arg(z). 



(A6) 
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